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ABSTRACT 



o 
o 
o 

Various analytic expressions have been proposed for the density profile of dark matter 
halos. We consider six of these expressions for which the density profile has a power-law 

■ fall-off p oc r~ 3 at large radii and a power-law cusp p oc r -7 (7 = 0, |, 1, |) at small 
radii. The phase-space distribution function for these models is calculated assuming 
spherical symmetry and either an isotropic velocity dispersion tensor or an anisotropic 
dispersion tensor of the type proposed by Osipkov and Merritt. The differential energy 

£T) '. distribution for these models is also derived. Several applications are discussed including 

the analysis of dark matter search experiments and the study of halo formation in a 
cosmological setting. Analytic fitting formulae for some of the models are provided. 

o " 

■ Subject headings: galaxies: halos — galaxies: kinematics and dynamics — dark matter 

6 

c/2 . 1. Introduction 

A fundamental problem in modern astrophysics is to determine the nature of relaxed systems 
^ ■ such as galaxies and dark matter halos. N-body methods are used extensively in this effort and 

allow one to follow explicitly the evolution of the phase-space distribution function (DF) /. A 
complementary approach seeks to obtain analytic steady-state models. Essential to this task is the 
Jeans theorem which states that the DF for any equilibrium system can be written in terms of the 
integrals of motion. In addition it is often possible to make an educated guess as to the form of the 
DF. Lynden-Bell (1967), Tremaine (1987) and Merritt, Tremaine & Johnstone (1989) for example, 
analyze model DFs that are designed to capture the physics of violent relaxation. By contrast, 
Henriksen &: Widrow (1999) propose a DF that is motivated by the spherical infall model wherein 
particle energies vary in a more orderly fashion. 

Unfortunately, neither the DFs considered by Merritt, Tremaine, & Johnstone (1989) nor the 
DF considered by Henriksen & Widrow (1999) lead to density profiles in agreement with those found 
in the simulations. The alternative is to begin with a desired expression for the density profile, 



1 widrow® astro .queensu.ca 



- 2 - 



together with simplying assumptions about the shape of the velocity ellipsoid, and construct the 
DF using standard techniques (e.g., Binney & Tremaine 1987, hereafter BT). Occasionally, one is 
lucky and finds an analytic DF-density profile pair, as with the Hernquist model (Hernquist 1990). 
Though for more general profiles such as the so-called "7-models" (see Carollo, de Zeeuw, and van 
der Marel (1995) and references therein) numerical integration is required, the problem is entirely 
tractible. 

The 7-models, which include the Hernquist model as a special case, have a density profile 
with a power-law cusp p oc r~ 7 at small radii and a p oc r~ 4 fall-off at large radii. While these 
models have proved to be useful in the study of elliptical galaxies, the dark matter halos found 
in cosmological simulations appear to have a power-law fall-off at large radii that is more gradual 
than r -4 . Navarro, Frenk, & White (1996) suggest that the density profiles of dark halos have a 
"universal" shape (the so-called NFW profile) of the form 

Po 

(r/a) (1 + r/af 

They obtain reasonably good fits using this expression for halos that range in mass from 3 x 10 11 M 
(dwarf galaxies) to 3 x 10 15 M (rich galaxy clusters). However, there is considerable debate over 
just what the profile is in the innermost regions of a halo. Kravtsov et al. (1998) find that the 
observed rotation curves of dwarf and low surface brightness galaxies can be fit by a profile with 
a shallower central cusp {p oc r~ 7 where 7 ~ 0.2) and their own N-body simulations support this 
conclusion. On the other hand Moore et al. (1998) have performed high resolution simulations of 
cluster-sized halos and find a central cusp that is steeper than r _1 (7 ~ 1.5). 



PNFW = ; — 7T2 W 



Kravtsov et al. (1998) advocate a general fitting formula of the form 

Cp 



(2) 



(r/a) 7 (1 + (r/af) (/3 " 7)/a 

7 controls the slope of the inner profile, (3 that of the outer profile, and a the sharpness of the 
transition. The normalization parameter C will be discussed below. The NFW profile corresponds 
to (a, (3, 7) = (1,3,1) while the 7-models correspond to (a, (3, 7) = (1,4,7) with the Hernquist 
profile appearing as the special case, 7 = 1. 

In this work, we derive semi-analytic DFs for a select subset of the models described by Eq.(2). 
Specifically, we focus on models with an r~ 3 density fall-off at large radii ((3 = 3). Furthermore, 
we consider only the six models for which the gravitational potential can be expressed in terms of 
elementary functions, namely a = 1; 7 = 0, \, 1, § and a = 2; 7 = 0, 1. 

In Section 2 the DFs for the six models assuming an isotropic velocity dispersion tensor arc 
calculated. DFs with anisotropic velocity dispersion of the Osipkov-Merritt type are considered 
in Section 3 and their existence and stability is discussed. In Section 4 the differential energy 
distribution various models is calculated. Several possible applications of our results are discussed 
in Section 5. Analytic fitting formulae for many of the DFs found in the text are provided in the 
Appendix. 
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2. Systems with Isotropic Dispersion Tensors 

1 /2 

For convenience we introduce the dimensionless variables R = r/a, V = vj {AirGp^a 2 ) , 
q = p/po and F = (AirG) 3 ^ 2 a 3 p^ 2 f. In addition we define the relative energy and relative 
potential (again in dimensionless form) to be respectively £ = — {E — $00) /AnGpoa 2 and \£ = 
- (<E> - $00) /AirGpoa 2 where E = \v 2 + <E> and <E> = <£(r) is the Newtonian potential with $(00) = 



- co- 



in general, the Newtonian potential calculated from the density profile (2) with [3 = 3 must 
be determined numerically (For the 7-models (a = 1; (3 = 4) an analytic form for the potential 
not only exists, but can be inverted to give r = r(&) in closed form.) However, the potential can 
be determined analytically for a = 1; 7 = 0, ^, 1, | (Models I-IV) and a = 2; 7 = 0, 1 (Models V, 
VI) and for convenience, we focus on these six cases. Expressions for the potentials are collected in 
Table 1. The normalization parameter C in Eq.(2) is chosen so that for \I>(0) = 1 with the limiting 
form asR—tOofty—tl — Ar 2 ~~< . Here A is a constant that depends on a and 7. For all models 
with an r~ 3 power-law fall-off at large radii, the asymptotic form as R — > 00 is ^ — > A' \n{R) / R. 

The distribution function for an equilibrium spherical system with an isotropic dispersion 
tensor depends only on the relative energy £ and can be calculated from the density profile and 
potential through an Abel transform (BT) 



F(£) 1 



d 2 g 1 / dp 

d^ 2 7£ =T ^ + 7^ v^y*= J 



V8tt 2 

In all of the models that we will consider, the second term on the right hand side is zero. 



(3) 



The integral that remains is evaluated numerically. The integrand diverges at one or both 
of the limits but this can be handled using standard techniques such as those found in Press et 
al. (1986). The DF is evaluated at values of £ equally spaced in ln£ for < £ < \ and equally 
spaced in In (1 — £ ) for \ < £ < 1. Accuracy is checked by calculating the density profile from the 
DF: 

q(R)=4V2tt I F {£) (*-£) 1/2 d£ (4) 
Jo 

and comparing with the exact expression. Typically, 2 x 10 4 integration points are required to 
guarantee 0.1% agreement over the range R = 1CT 6 — 10 6 . 

In Figure 1, we compare DFs for the Hernquist and NFW profiles. The DFs in the regime 
(1 — £) -C 1 are nearly identical (up to a normalization constant). This is to be expected since 
in this regime the systems are dominated by particles at small radii and the Hernquist and NFW 
profiles each have an r _1 cusp. While the DFs diverge in this limit the mass at small radii is finite, 
as we will see in Section 5. 

The DFs for both the Hernquist and NFW models decrease as £ — > 0. However, the decrease 
is slower in the NFW model, a reflection of the fact that the halo in this model is more extended. 
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It is straightforward to determine the functional form of the DF in this limit. For the NFW model 

d 2 g 1 * 



d^ 2 R (lni?) 2 (-In 

and we find 



(5) 



^3/2 

as compared with the Hernquist model where F oc £ 5 / 2 . 



F (£) « ; ; „ n3 (6) 



The results for our series of a = 1 models (I- IV) are plotted in Figure 2. All of the DFs in this 
series have a limiting form as £ —>■ given, up to a constant, by Eq.(6). The difference in the index 
for the power-law cusp at small radii is reflected in the behavior of the DFs as £ — ► 1 with the 
trend that a steeper inner cusp corresponds to a stronger divergence in this limit. Once again, we 
can determine a limiting form for the DF, this time as £ — ► 1. For 7 = (77, 1, |), d 2 p/d^> 2 oc i? 7-4 
and we find 

F ( £ ) * f° — Vh (7) 

OC (1 - f )-(6-7)/(4-2 7 ) (8) 

where i?£ is defined by the relation ^(Rg) = £. In writing this expression, we use the fact 
that the integral is dominated by the region in R ~ R £ oc (1 - £) 1/{2 ~ 7) . For the NFW profile, 
F oc (1 — £)~ 5 ^ 2 which is the same as is found in the Hernquist model (Hernquist 1990). The case 
7 = is handled separately: We find d 2 p/d^ 2 oc R~ s and F oc (1 — £ 

Kravtsov et al. (1998) suggest that the plausible value of the parameter a = 2 corresponding to 
a sharper transition between inner and outer regions of the halo. In Figure 3 we compare the DFs 
for the two models (a, (3, 7) = (1, 3, 0) and (2, 3, 0). The most striking difference occurs as £ — > 1 
(small radii) where the DF for the a = 2 model approaches a constant. This can be understood by 
noting that as R — > 0, d 2 g/d^ 2 is finite for a = 2; 7 = but diverges for a = 1. 

The distinction between the DFs for (a, /3, 7) = (1, 3, 1) and (2, 3, 1) is more subtle. The 
asymptotic forms as £ — > 1 and £ — > are the same and so the difference between the models arises 
solely in the transition region. In Figure 4 we plot the ratio of the DFs for the two models. If we 
normalize the models at £ = 1 (R = 0) then as R — > 00 the DF for a = 2 will exceed that for a = 1 
by a factor ~ 6. 

Physical DFs (F > for < £ < 1) exist for all of the isotropic models considered here. 
Moreover, both F and d 2 g/d^ 2 are monotonically increasing functions of £. This, by Antonov's 
second and third laws, is sufficient to guarantee that these models are stable against both radial 
and nonradial perturbations (BT). 
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3. DFs with Anisotropic Velocity Dispersion 



The DF of any steady-state system that is spherically symmetric in both velocity and con- 
figuration space can be expressed as a function f(E,L) where L is the magnitude of the angular 
momentum vector (BT). In general, the velocity-dispersion tensors for these models are anisotropic: 
the velocity dispersions in the two tangential directions are equal but in general different from the 
velocity dispersion in the radial direction. We investigate a special class of these models in which 
the DF assumes the form (Osipkov 1979; Merritt 1985a,b) 



where £ = L/ (4irGp a 4 ) 1 and Q = £ - C 2 /2R 2 a . In addition, the condition F = for Q < 
is imposed. R a is often called the anisotropy radius: Outside R a the velocity dispersion is peaked 
toward radial orbits while inside R a the dispersion is nearly isotropic. As discussed in Carollo, 
de Zeeuw, and van der Marel (1995), physical models do not exist for R a < R ajia - m where -R a ,min 
depends on the shape of the potential and must be determined numerically. For values of R a below 
this threshold, the DF must become negative for a range of values in Q in order to compensate for 
the excess population of radial orbits needed to produce the halo at large radii. The problem is more 
severe for the systems considered here since the density profile at large radii varies as r~ 3 rather 
than r -4 . In particular, we find that R a ,mm — 0.36 and 0.75 for (a, (3, 7) = (1, 3, 0) and (1, 3, 0) 
respectively. By comparison, Carollo, de Zeeuw, and van der Marel (1995) find i? a ,min — 0.20 and 
0.45 for (a, /3, 7) = (1, 4, 0) and (1, 4, 0) respectively (cf. their Figure 1). 

The DF for the Osipkov-Merritt models is found by replacing the g in Eq. (3) with the auxiliary 
density qq = (l + R 2 /R 2 ) g. The results for a sequence of Osipkov-Merritt NFW models is shown 
in Figure 5. The pathological nature of the DF as R a — ► i? a ,min is evident in the dotted curve 
(R a = R a ,min + e where e ~ 10" 5 ). 

Notice that for these models the DF at small £ is a decreasing function of £ (F oc £~ l l 2 ). 
Stability analysis against radial perturbations by application of Antonov's second law is therefore 
inconclusive (BT). While models with R a — ► 00 are stable against both radial and non-radial 
perturbations those with R a — ► are almost certainly unstable to radial perturbations (Merritt 
1985b). Numerical experiments are therefore required to determine the exact region of stability. 



When comparing results of N-body simulations with those from analytic models, it is natural 
to use the differential energy distribution dM/d£ (BT). For models with an isotropic velocity 
dispersion tensor, dM/d£ is simply the product of the DF with the density of states: 



F(£,C)=F(Q) 



(9) 



4. Differential Energy Distribution 




(10) 
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where 

[■Re 

G{£) = 16\/2vr 2 / - £) 1/2 R 2 dR (11) 
Jo 

M is a dimensionless mass variable scaled by poa 3 . Note that as £ — ► 0, G(£) oc £~ 5 / 2 (up to 
logarithmic corrections) for both the models considered here and the Hernquist model. 

Figure 6 is a plot of dM/d£ for models I-V and for the Hernquist model. In the limit £ — > 0, 
dM/d£ for the Hernquist model approaches a constant reflecting the fact that the total mass is 
finite. Conversely, for the models considered here, dM/d£ diverges as £ -1 , a symptom of the 
logarithmic mass divergence at large radii. 

It is a somewhat more difficult exercise to calculate the differential energy distribution for 
Osipkov-Merritt models. Let us begin, as is done in the derivation of Eq.(10), with the expression 
for the total mass within a radius R (BT): 

r-R /•*(#') /•t/2 

M(R) = 16ir 2 R /2 dR' vd£ F(£, RV sin 77) sinjydr? (12) 

Jo Jo Jo 

where 77 is the polar angle in velocity space: V r = V cos r/. We first replace the integration variable 77 
with Q using the relation Q = £ — \ (R/R a ) 2 V 2 sin 2 rj and then interchange the R and £ integrations 
to find 1 r £ 

M(R) = 8V2ir 2 R a C d£ [ £ R'dR' I F (Q) d Q (13) 

Jo Jo Je-e (Q- {£-£)) ' 

where £ = \ (R' /R a ) 2 V 2 . Thus, the differential energy distribution for the Osipkov-Merritt models 
can be evaluated by performing the double integral 

™ = f""" km f FiQ)dQ m <i4) 

<l£ Jo Je-e (Q - (f - f)) 1/2 

where R m i n = min(R£, R). It is straightforward to show that in the limit R a — ► 00, £ — > and 
Eq.(14) reduces to (10). 

The results for the four NFW models considered in Section 3 are shown in Figure 7. Note that 
the divergence as £ — ► is even more severe (dM/d£ oc £~ 2 ). 



5. Applications 

In this section three possible applications of the results presented above are briefly discussed. 



5.1. Initial Conditions for N-body Experiments 



The dynamic range in N-body simulations has now reached the level where it is possible to 
study individual halos in exquisite detail while maintaining an accurate representation of the large- 
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scale tidal fields that shape them. Nevertheless, it is often useful to conduct highly controlled, 
albeit artificial, numerical experiments in order to gain a better understanding of the processes at 
work during halo formation. For example, numerous groups have used numerical simulations to 
study the accretion of satellites with a parent halo as well as the merger of two halos of comparable 
mass. Another set of experiments follows the evolution of a gas and dark matter system inorder 
to better understand the process of disk formation. Typically, very simple DFs (e.g., truncated 
isothermal sphere, Plummer model) are used to set up the initial conditions for these experiments 
and the choice is often made out of convenience, i.e., the availability of analytic DFs, rather than 
an expectation of what might arise in a cosmological setting. 

The DFs found in this work allow one to set up a variety of halo models with various forms for 
the density profile and velocity distribution tensor. For the models with isotropic velocity dispersion 
tensors, the probability for a particle to have energy £ and radius R is 

P(£,R) oc R 2 (V - £) 1/2 F (£) (15) 

The rejection method described, for example, in Press et al. (1986), provides a simple technique for 
selecting particles from this distribution. Once £ and R for a given particle are known, the speed 
V can be determined immediately. One then chooses, at random, four angles (two in configuration 
space and two in velocity space) to yield the six phase space coordinates of the particle. This 
procedure can be extended easily to the Osipkov-Merritt models. 

5.2. Interpretation of Results from Dark Matter Search Experiments 

The announcements in 1993 by the MACHO (Alcock et al. 1993) and EROS (Aubourg et 
al. 1993) collaborations of candidate microlensing events toward the LMC have highlighted the need 
for self-consistent model DFs of the Galaxy's halo. Interpreting the results from these experiments 
requires a comparison of the observed and predicted event rates. The latter has now been calculated 
for a set of halo models too numerous to list here. These models are generally constructed in one 
of two ways. One can begin with an ansatz for the mass distribution (e.g., triaxial spheroid with a 
prescribed density law) and assume a simple form for the velocity dispersion tensor (e.g., isotropic 
and Maxwellian). Of course, models constructed in this manner do not, in general, correspond to 
true equilibrium systems, i.e., solutions of the time-independent collisionless Boltzmann equation. 
A second approach employs exact, analytic model DFs. Here, the Jeans theorem is invaluable 
since it allows one to construct model DFs simply by taking functions of the integrals of motion. 
The so-called power-law models in which the DF is constructed from powers of the energy and the 
angular momentum vector L z (Evans 1993) have been used in this way (Evans Sz Jijina 1994). 

The DFs presented in this work provide an alternative set of models suitable for the analysis and 
interpretation of results from dark matter search experiments. The density profiles include many 
of the popular forms found in the literature. In addition, the Osipkov-Merritt ansatz for velocity 
space anisotropy is consistent with what one expects for a realistic halo model, i.e., primarily radial 
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orbits in the outermost regions of the halo where particles have only recently separated from the 
Hubble flow. 

5.3. Formation of Dark Matter Halos 

The hierarchical clustering hypothesis provides a compelling picture for the formation of struc- 
ture in a Universe dominated by cold, dissipationless dark matter. This scenario naturally lends 
itself to the development of phenomenological models for the growth of dark matter halos (e.g., the 
Press-Schechter formalism (Press and Schecter 1974)). Nevertheless, hierarchical clustering says 
little about the internal structure of these systems. For this, one must understand the process 
by which a collisionless system relaxes to an equilibrium or quasi-equilibrium state. To this end, 
two opposing pictures have emerged. The first suggests that violent relaxation is the dominant 
process at work, the conjecture being that particles will transfer energy to one another as they 
move through the rapidly varying potential of the collapsing system. In so doing, the particles 
lose all memory of their initial state. The alternative picture is based on the spherical infall model 
(Gunn and Gott 1972; Fillmore and Goldreich 1984; Bertschinger 1985) in which a system relaxes 
from the inside out. Particles near the peak of the initial density perturbation collapse first and 
constitute the most tightly bound regions of the final system. Likewise, particles in the wings of the 
initial density perturbation collapse later on and form the system's outer halo. Here, one expects a 
direct correspondence between initial conditions and the final state of the system and in particular, 
a tight correlation between initial and final energy and angular momentum. 

Not surprisingly, the dark matter halos found in cosmological simulations appear to follow 
an intermediate path in reaching a relaxed or virialized state. Quinn and Zurek (1988) track the 
binding energy and angular momentum of selected particles during the simulated formation of a 
dark matter halo. They conclude that though energy and angular momentum are not conserved, 
the collapse is more orderly than one would predict assuming violent relaxation. Indeed there seems 
to be a tight correlation between initial and final energy and angular momentum. The suggestion 
then is that there is a simple relation between the differential energy and angular momentum 
distributions of a relaxed system and the initial perturbation that gave rise to it. The DFs presented 
in this work may provide a further link between the dM/d£ or d 2 M/d£dL and the density profile. 

6. Conclusion 

The density profile and associated gravitational potential provide a popular way of character- 
izing dark matter halos. To be sure, the halos found in simulations are varied in shape and rich in 
substructure. However, it is essentially the density profile that determines the contribution by the 
halo to the observed rotation curve. Moreover, the evolution of the gas component (e.g., in forming 
a disk galaxy) depends sensitively on the shape and depth of the potential. 
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In this paper, we have presented DFs that reproduce six different density profiles. All of the 
models assume spherical symmetry and a power-law fall-off at large radii p oc r~ 3 but differ in the 
slope of the inner power-law cusp and the sharpness of the transition in going from inner to outer 
regions of the halo. 

We have considered models with velocity-space anisotropy of the type proposed by Osipkov and 
Merritt. The key feature and also limitation of these models is that the distribution of orbits varies 
from one that is nearly isotropic in the inner regions of the halo to one that is nearly radial in the 
outer regions. In addition, since the Osipkov-Merritt models provide a family of DFs, parametrized 
by R a , that reproduce the same density profile, they may be combined to yield a fairly general class 
of models (Merritt 1985a). 

Though the DFs presented in this work were derived numerically the analytic fitting formulae 
provided in the Appendix should enable researchers to use them with the same ease as they would 
find with closed form DFs such as the Hernquist model. 

I would like to thank the Canadian Institute for Theoretical Astrophysics for their hospitality 
during a sabbatical visit. I also acknowledge R. Henriksen for useful discussions. This work is 
supported by the Natural Sciences and Engineering Research Council of Canada 
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A. Appendix 

In this appendix, we provide analytic fitting formulae for some of the DFs presented in the 
text. The limiting forms for the isotropic models, as discussed in the text, are: 

m J f /2( -'f r3 f r £ r° (ad 

V ' \ (1 - £y x for £ -> 1 v ' 

where A = (l, ^, §, §, 0, §) for Models I- VI respectively. Noting that -ln£ ~ (1 - 5) for l-£ <C 
1, we propose the following fitting formulae for the isotropic DFs: 



F(£) = F £ 3 ^ (1 - £)~ x (^|^) 9 e P (A2) 



where the polynomial P = YliPi^ 1 IS introduced to improve the fit. In general, 3-5 terms in P are 
required to achieve a reasonable fit. The parameters are given in Table 2. Notice that q ~ — 3 as 
is expected from Eq.(Al). 

For the anisotropic models discussed in the text (a = 1; 7 = 1), we use the fitting formula 

F(Q) = F Q-^ (1 - Q)' x (^p) 9 e P (A3) 

where P is now a polynomial in Q. The parameters are given in Table 3. No attempt is made to 
fit the DF for R a ~ i? a ,min- 



Table 1. Newtonian Potentials for Selected Density Profiles 



Model a 7 C * 



I 


I 





2 


|ln(I+ J R)-(I+ J R)- 1 


II 


I 


1 

2 


3 
2 


1 (i?+^)i/2 + 2R lnS 


III 


I 


1 


1 




IV 


1 


3 
2 


1 
2 




V 


2 





1 


^sinn" 1 ^) 


VI 


2 


I 


2 

7T 


l_l tan -i (i?) + J_ ln ( 1 + jR 2) 



Note. — S = 1 + 2R + 2 (R + i? 2 ) 1 



Table 2. Fitting Formula Parameters for Isotropic Models 



Model 


Fo 




q 


Pi 


P2 


P3 


Pa 


P5 


P% 


I 


2.0460 x 10- 


-2 


-2.7145 


-1.0215 


23.766 


-98.330 


194.50 


-180.01 


63.296 


II 


3.6478 x 10" 


-2 


-2.7092 


0.8670 


-10.035 


65.895 


-166.31 


179.21 


-70.007 


III 


9.1968 x 10" 


-2 


-2.7419 


0.3620 


-0.5639 


-0.0859 


-0.4912 






IV 


4.8598 x 10" 


-1 


-2.8216 


0.3526 


-5.1990 


3.5461 


-0.8840 






V 


5.8807 x 10" 


-2 


-2.6312 


-3.7147 


41.045 


-132.20 


216.90 


-170.23 


51.606 


VI 


1.4696 x 10" 


-1 


-2.6210 


-3.6125 


23.172 


-78.104 


135.80 


-123.11 


43.705 



Table 3. Fitting Formula Parameters for Anisotropic Models 



Ra 


F 






Pi 


P2 


P3 


P4 


P5 


P6 


0.6 


1.0885 x 10" 


-1 


-1.0468 


-1.6805 


18.360 


-151.72 


336.71 


-288.09 


85.472 


1 


3.8287 x 10" 


-2 


-1.0389 


0.3497 


-12.253 


-9.1225 


101.15 


-127.43 


47.401 


3 


4.2486 x 10" 


-3 


-1.0385 


0.7577 


-25.283 


149.27 


-282.53 


229.13 


-69.048 


10 


3.8951 x 10" 


-4 


-1.0447 


-2.2679 


79.474 


-237.74 


329.07 


-223.43 


59.581 
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Fig. 1. — DFs as a function of relative energy 6 for the NFW profile (solid curve) and Hernquist 
profile (dashed curve). 
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Fig. 2.— DFs for Models I-IV (a = 1; j3 = 3): Model I, 7 = (dotted curve); Model II, 7 = \ 
(long-dashed curve); Model III, 7 = 1, the NFW profile (solid curve); Model IV, 7 = § (dot-dashed 
curve) . 
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Fig. 3. — Comparison of the DFs for Models I and V (7 = 0; (3 = 3): a = 1 (dotted curve) and 
a = 2 (dot-long dashed curve). 
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Fig. 4.— DF for Model VI divided by the DF for Model II 
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Fig. 5. — DF for the NFW profile assuming a velocity dispersion tensor of the type proposed 
by Osipkov and Merritt. The different curves correspond to different values of the anisotropy 
parameter R a : R a = 0.35548 (dotted curve); R a = 0.6 (dashed curve);i? a = 1 (long-dashed curve); 
R a = 3 (dot-dashed curve); R a = 10 (dot-long dashed curve). The DF assuming an isotropic 
dispersion tensor R a — > oo (Model III) is given by the solid curve. 
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Fig. 6. — Differential energy distribution, dM/d£ as a function of energy for Models I-V and the 
Hernquist model. Curves are labeled as in Figures 1 and 2 
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Fig. 7. — dM/d£ for the anisotropic models. Curves are labeled as in Figure 5 




